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Abstract 

Lack of Lorentz invariance of QCD at finite quark chemical potential in gen- 
eral implies the need of Lorentz non-invariant condensates for the self-consistent 
description of the color-superconducting ground state. Moreover, the spontaneous 
breakdown of color SU (3) in this state naturally leads to the existence of 5f7(3) 
non-invariant non-superconducting expectation values. We illustrate these observa- 
tions by analyzing the properties of an effective 2-flavor Nambu-Jona-Lasinio type 
Lagrangian and discuss the possibility of color-superconducting states with effec- 
tively gapless fermionic excitations. It turns out that the effect of condensates so 
far neglected can yield new interesting phenomena. 



1 Introduction 



According to current wisdom the deconfined 2-flavor QCD matter at moderate baryon 
densities and low temperatures behaves as a color superconductor Standard charac- 
teristics of its BCS-type ground state is the ground-state expectation value p|, ^ 

5i = (V'^C75r2A2^) , (1) 

which corresponds to a scalar diquark condensate in a color anti-triplet state. Here the 
superscript T denotes a transposition and C the matrix of charge conjugation. T2 is a 
Pauli matrix acting in flavor space, A2 a Gell-Mann matrix acting in color space. We 
have taken the freedom to rotate the unpaired quark color into the 3-direction, which can 
always be done. 

Because of the empirical fact that the (approximate) chiral SU(2) symmetry of the 
QCD Lagrangian is not respected by the QCD vacuum, it is natural to ask whether the 
quark (-antiquark) condensate 

= (V^V^), (2) 

persists also in the ground state of QCD matter at finite baryon density. For model 
Lagrangians with exact chiral symmetry one usually finds a first-order phase transition 
from a low-density phase (or vacuum) with 7^ and (5i = to a high-density phase 
with 5i 7^ and = 0. This is different if there is a small quark mass m which explicitly 
breaks chiral symmetry. In this case cannot exactly vanish above the phase transition 
and coexists with the diquark condensate. In fact, just above the phase transition the 
gaps related to the two condensates can be of similar magnitude Q. 

At finite density the existence of Lorentz non-invariant expectation values becomes 
possible. The most obvious example is of course the density itself, 

P = (^tV), (3) 

which transforms like the time component of a 4-vector. Consequently, if the Pauli prin- 
ciple and the form of the interaction permit, the ground state of a 2-fiavor color supercon- 
ductor is characterized by more condensates than merely the Lorentz-invariant ones, (|I|) 
and d^). This simple observation was made already in the pioneering work of Bailin and 
Love p. For instance, there could be another diquark condensate of the form [|, |, 

S2 = (V^^C7°75T2A2^) , (4) 

which also transforms like the time component of a 4-vector. The ground state does not 
even have to be isotropic, but rotational invariance could be spontaneously broken, like 
in a ferro magnet. Even the breakdown of translational invariance due to the formation of 
crystalline phases is a possible scenario P, P). 

Up to now, we have discussed only space-time symmetries. We should recall, however, 
that in the presence of the diquark condensates (P or (^) also color SU{3) is broken, 
since only the first two colors ("red" and "green") participate in the condensate, while 
the third one ("blue") does not. Therefore, there is no reason to assume that all other 
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condensates are co\ot-SU{3) invariant in this state. For instance, we should expect that 
the contributions of red and blue quarks, (pr and (pb, to the quark condensate could be 
different, thus giving rise to a non- vanishing expectation value 

08 = (^Ag^) = (5) 

(Note that 6i and 62 leave a color SU{2) subgroup invariant and therefore all green 
quantities are identical to the red ones). Similarly, the densities of red and blue quarks 
will in general not be the same, i.e., in addition to the total number density p = 2pr + pb 
there could be a non-vanishing expectation value 

2 

P8 = (^7°A8^) = -^{Pr-Pb)- (6) 

Since these color-symmetry breaking expectation values, induced by the presence of color- 
symmetry breaking diquark condensates, could in turn influence the properties of the 
diquark condensates, in principle, all condensates should be studied in a self-consistent 
way. 

Detailed understanding of the color-superconducting state which in Nature might ex- 
ist in the interiors of neutron stars [110] requires, however, a detailed knowledge of the 



effective quark-quark interactions close to their Fermi surface [|ri|] . Due to the lack of ex- 
perimental data or information from the lattice such a knowledge is missing at present. It 
is nevertheless possible to analyze in detail the effective low-energy quantum field theory 
of the deconfined low (or zero) temperature quark matter within models which respect 
all relevant symmetries of the corresponding QCD Lagrangian. Lack of Lorentz invari- 
ance of QCD at finite p implies that Cejf itself should respect only an 0(3) rotational 
symmetry. Since there are no thermal gluons at T = it is justified and customary to 
analyze the ground-state properties of deconfined quark matter by virtue of Cejf having 
only a global color SU{3) symmetry. Response of gluons to the quark condensates car- 
rying color (Meissner effect) is then studied |13| perturbatively, i.e., as if the color 
gauge symmetry is spontaneously broken (for more detailed discussion see [Q). In this 
article we illustrate the need of Lorentz- and color-5'[/(3) non- invariant condensates for 
the self-consistent description of the color-superconducting ground state. For simplicity 
we leave the rotational symmetry unbroken, i.e., we restrict ourselves to the isotropic case. 

The remainder of this article is organized as follows. In Sec. ^ we will derive the 
thermodynamic potential and a coupled set of gap equations for a general NJL-type model, 
taking into account a possible condensation in the channels (|l]) to (H). The resulting 
dispersion laws are discussed in more detail in Sec. ^. Here special emphasis is put 
on possible gapless color superconductors. In Sec. § we discuss our numerical results. 
Conclusions are drawn in Sec. |^. 
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2 General Formalism 

In order to illustrate the generic properties of the superconducting ground state let us -as 
an example- consider a model defined by the Lagrangian density 

^eff = i'{i^-m)i) + Cqq + Cgq (7) 

describing a quark of mass m which interacts with other quarks and antiquarks via NJL- 
type 4-point interactions of the form 



-qq 

and 



A.- = #(^^)' + 9?H^^a^? + gfH^i^f + gfHi'i'Ki^f + - (8) 



Cqq = /il (V'n5^2AA^'')(^''^75T-2AAV') + ^2 7 75^2Aa V''') (V''' 7 75^2 ^) + ••• (9) 

Here the dots indicate possible other channels, not related to the expectation values (|l|) 
to d^). In Eq. (y) we used the notation ip'^ = Cip^ , while \a denotes the antisymmetric 
Gell-Mann matrices A2, A5, and A7. All color indices are understood to be summed over. 
In the following it is convenient to formally double the degrees of freedom by defining 

= 4. flS) . (10) 



^2 V^'(^) 

To obtain the mean-field thermodynamic potential at temperature T and chemical poten- 
tial /i we linearize £e// in the vicinity of the expectation values Eqs. (|I|) to and apply 
Matsubara formalism. The resulting thermodynamic potential per volume reads 

^{T,^,) = -TY, I ^ ^Tr \n(^^S-\zu;^,p)) + V, (11) 

with Un being fermionic Matsubara frequencies. S~^{p) is the inverse propagator of the 
q-fields at 4-momentum p. It is given by 



^- Mq- MgAg + /io7° + m^Ag Ai75r2A2 + A27°75i"2A2 
-A*75r2A2 + A;7075r2A2 ^ - Mo - MgAg - /io7° - f^sl'^^s 



(12) 



Here we introduced the effective quark masses 

Mo = m-2gi'U , Mg = -2gi'Us , (13) 
the effective chemical potentials 

and the diquark gaps 

Ai = -2hi6i , A2 = 2/12^2 . (15) 
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These quantities also enter into the last term of Eq. (|TT]), which is defined as 

^ (Mo-m)^ , Mt , (/io-/i)' , JfL : ^ : ^ rifi^ 

For later convenience, but also for the interpretation of the results, it is useful to perform 
linear combinations to get red and blue quantities, e.g. red and blue constituent quark 
masses = Mq + ^ Mg and = Mq - ^ Mg. We then find 

= m - ^(6^f + - ^(3# - 2gf)c^, , 

M,=m- ^(6# - 4^f - ^(3# + 4^f , 



= /i + ^(6^?f - 4^?(«))p, + ^(3^?f + 4^?(«))p, . (17) 

For two fiavors and three colors the inverse propagator, Eq. ([T2|) , is a 48 x 48 matrix, 
and the trace in Eq. (^) has to be taken in this 48-dimensional space. After performing 
the Matsubara sum we obtain: 

0(T, /^) = -4 / { 2 [^±±^ + T HI + e---/-) + T ln(l + e"-/-)) 

+ ( eb + T ln(l + e-^+/^) + T ln(l + e-^-/^))} 
+ V, (18) 

where physically irrelevant constant terms have been suppressed. The dispersion laws 
which enter into this expression are given by 



e± = efcipb = ./p2 + M2±/ib (19) 



and 



with 



E± = v'p^ + M2 + /i2+ |Ai|2 + |A2|2±2s, (20) 



s = V(p2 + |A2P)p2 + t2^ t = Mri^r - Re{AiA*) . (21) 

So far, the thermodynamic potential depends on our choice of the expectation values 
(|l|) to (^, which determine the effective masses, effective chemical potentials and diquark 
gaps as indicated above. On the other hand, in a thermodynamically consistent treatment 
the condensates should follow from the thermodynamic potential by taking the appropri- 
ate derivatives. The self-consistent solutions are given by the stationary points of the 
potential, 

6n 6n 6n 6n 6n 6n 

TTT = TTT- = 1— = 1— = 1^ = 1^ = 0. 22 

6Mq 6Ms 6ho 6fis SAi 6A2 
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If there is more than one stationary point, the stable solution is the one which corresponds 
to the lowest value of Q{T,^). 

Eqs. (0) and ( ^2]) lead to the following expressions for the various expectation values: 



= -4 



% = 



Pr = 



+ (1 - 2n{E. 

J (27r)3 2s I ^ ^ ^ 

+ (1 - 2n{E_ 



pb 
5 



1 — 



4 
-4 



(n(e_)-n(e+)) 

d?p 1 
(27r)3 I 



(i^p 1 



(27r)3 s 



(1 - 2n(E+ 
+ (1 - 2n{E^ 
(1 - 2n(^+ 
+ (1 - 2n{E^ 



E+ 

— [MrS - flrt] I , 
^ [firis + f) + Mrt] 

E+ 

-^[^,{s-f)-Mrt] } 



-^[Ais + Ast] }, 

^[A2(s+p-^)-Ait] 

E+ 

-^[A2(s-/) + Ait] 



(23) 



where n(E) = l/(e-^/^ + 1) is a Fermi function. Together with Eqs. (|15|) and (|1^ these 
equations form a set of six coupled gap equations for M^, M;,, /x^, /i^, Ai and A2. The 
expressions for the blue expectation values 0b and formally look like the corresponding 
formulae for free particles. However, they depend on the effective quantities Mf, and 
which are also related to red quantities via Eq. ([T^). Despite of these interdependencies, 
the masses of red and blue quarks, and also their densities will in general be different, as 
anticipated above. 

Another interesting observation is that in general 5i and 82 or, equivalently, Ai and 
A2 cannot vanish separately. This means that the familiar scalar diquark condensate 5i 
will in general be accompanied by an induced non- vanishing expectation value 82- Q On 
the other hand, there is always a solution with Ai = A2 = 0. In this case the expressions 
for 0r and get the same structure as the analogous expressions for the blue quarks. If 
Ai and A2 do not vanish, they are determined by the gap equation only up to a common 
phase. In most cases this phase can be chosen such that Ai and A2 are real. 



*This point has already been discussed in Ref. ||7|. The authors, however, argued that because of the 
small value of §2 as compared with 5i , this condensate could for simplicity be neglected although strictly 
speaking in this case the system of gap equations cannot be closed. 



6 



3 Dispersion laws and gapless color-superconducting 
states 



In physical terms the derived mean-field thermodynamic potential, Eq. ([Tq), describes the 
thermodynamics of a mixture of noninteracting fermionic excitations of two types: (A) 
quark excitations with the dispersion law e± and (B) Bogoliubov-Valatin (BV) quasiquark 
excitations with the dispersion law Straightforward interpretation is to say that these 
dispersion laws should be used for calculating (in principle) measurable quantities, e.g. 
the specific heat of the system. The parameters entering the dispersion laws, i.e., Mb, 
Mr, fib, /^r, l^il; l^2|, and cos(v9i — (^2), i^j = l^il exp (2(y9j), j = 1,2) are fixed as 
the solutions of the coupled nonlinear integral equations (p3D in terms of m, /x, and the 
dimensional coupling constant (s) of a given four-fermion interaction. Due to the fact that 
the integrals in Eq. (^3]) have to be regularized there is one more dimensionful parameter 
upon which the above quantities depend. This can be a cutoff, a parameter in a form- 



factor [g], the instanton size [0, or a thickness of a layer around the Fermi surface in 
which the interaction is assumed to be different from zero 0]. 

Component (A) of the system behaves like a normal relativistic Fermi gas, character- 
ized by a linear dependence of the low-temperature specific heat on T. The actual be- 
havior of the physically interesting BV component of the system characterized by (p) 
depends strongly upon the details of the interaction. For instance, the instanton mediated 
Lagrangians employed in Refs. [0] and P, IT^ have /12 = and therefore A2 = 0. In this 



case Eq. (g^) reduces to the "classic" result 

E- (p) = ^(v/j^ + M2-/i,)2 + |Ai|2 . (24) 

There are no nodes in this function, and this component exhibits a superconducting 
behavior including the exponential low-T specific heat characteristic of ordinary super- 
conductors. 

In the general case the form of E-{p) is given by Eq. (pO|). Clearly, without Lorentz 
invariance the dependence of the energy of a particle-like excitation on momentum is 
restricted only by positivity. It is not surprising, however, that E_{p) can be parametrized 

as 



E^P) = ^i^P" + Miff - Mf? + |Ae//P , (25) 

provided we identify 

Mf = y ^^r + \^2\, Kff = 71 , (26) 



and 



Mr' + \A,\' = M^jf+\Aeff\\ (27) 



'I'Strictly speaking, the quarks of component (A) are also quasiparticles as long as Mi, is different from 
the bare mass m. 
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It is interesting to notice that |Ae//| can vanish, i.e., E-{p) can have nodes at 
/igjj — Mlj^. This is the case if 



If /ir 7^ this relation can be used to ehminate Ai in Eq. (^): 



/^e// = (1 + ^) , ^e// = M,2(l + ^), (29) 



and hence 



fnode = if^r " M^) (1 + ^) . (30) 

This means, /i^ must be greater or equal to and it immediately follows from Eq. (|28|) 
that a gapless color-superconducting solution is only possible if IA2I > |Ai|. 

In the vicinity of the node the BV quasiparticle takes the form of a non-relativistic 
fermion, 

E_{p)^^-f,*, (31) 
2m* 

where m* = fi^ff and /i* = |(/ie// — 7777)' Despite the superconducting condensates the 
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specific heat of such a system is linear in T. A similar phenomenon was found in Ref. 
in the color-flavor locked (CFL) phase, though its origin is not the same as here. 

In the above discussion it is tacitly assumed that there is an interaction which yields 
solutions of the coupled gap equations for which Eq. (^) holds. A particularly simple way 
to fulfill this condition is to assume that M,. = Ai = 0, which could be realized by taking 
m = gf''' = hi = and only /i2 7^ 0. If we regularize the divergent integrals by a sharp 
3-momentum cut-off A and restrict ourselves to T = , the thermodynamic potential of 
this schematic model is readily calculated: 



n 



schem(T = 0,/i;A2) = -^(2(^2 + |A2|')' + /i') + (32) 



Here we dropped an irrelevant constant —^A^. Obviously this function is not bounded 
from below, but only the self-consistent solutions of the gap equation, i.e. 6Qschem/ SA2 = 
0, are physically meaningful. There is always a trivial solution with A2 = 0. For < 
^2 < 1^ there are also nontrivial solutions with IA2P = — /i^. However, whenever 
these nontrivial solutions exist, they correspond to maxima of ^Ischem, while at same time 
the trivial solution is a local minimum with a lower value of ^Ischem- This means, the 
nontrivial gapless solution is unstable. Although we have shown this only for our very 
simple schematic model, it might be a rather general feature. In fact, a similar observation 
was made by the authors of Ref. [|1^ for the CFL phase. We will come back to this point 
in the end of the next section. 
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4 Numerical results 



In the schematic example discussed in the end of the previous section the problem was 
reduced to a single condensate. In this section we want to present the results of a numerical 
study of the full coupled set of gap equations derived in Sec. 0. As an example we consider 
an interaction with the structure of a "heavy-gluon exchange" , 

Chge = i^{i^-m)tp - gE{ip1^\ai^f + Qm {'^P^Ki^Y , (33) 

although there is no reason why the effective interaction at moderate densities should 
have this particular form. As discussed in the introduction, the effective Lagrangian at 
finite densities does not need to be Lorentz invariant. We underline this possibility by 
explicitly allowing for different "electric" and "magnetic" coupling constants, qe and qm- 
The effective quark- antiquark interaction Cqq and the effective quark-quark interaction 
Cqq as given in Eqs. (||) and @ can be derived from Chge by performing the appropriate 
Fierz transformations. The resulting coupling constants are 

2 1 1 

gf^ = g + 3^m) , gf^ = -— {gE + 3^m) , hi = -{gE + Sqm) , 

2 11 

= n^9E- SgM) , gi^^ = -gE - — {gE- Sgu) , h2 = -{qe- SgAi) ■ 
y 24 D 

(34) 



4.1 Lorentz invariant interaction 

We begin our discussion with the case of a Lorentz-invariant interaction, 

g = gE = ghi ■ (35) 

Although there is no reason for this assumption, it allows for a better comparison with 
other calculations in the literature, which often start from an interaction of this form. Of 
course, having a Lorentz-invariant interaction does not mean that there are only Lorentz- 
invariant condensates, since Lorentz-invariance is still broken by the chemical potential. 

If we insert Eq. (^) into Eq. ( P^ we find that for g > the interaction is attractive 
in the scalar quark-antiquark channel and in the scalar diquark channel and repulsive in 
all other channels of interest. Of course, non- vanishing expectation values in the repulsive 
channels do not develop spontaneously, but only as a result of an external source, like the 
chemical potential, or induced by non- vanishing expectation values in attractive channels. 
In fact, the solutions of the gap equations correspond to maxima of the thermodynamic 
potential with respect to variations in the repulsive channels, whereas they can be maxima 
or minima with respect to variations in the attractive channels. 

In our numerical calculations we restrict ourselves again to T = and take a sharp 
3-momentum cut-off A to regularize the integrals. Although a quark model of the present 
type should be used only in the deconfined phase, we follow the general custom and fix 
the parameters such that they yield "reasonable" vacuum properties. In the following 
we take A = 600 MeV, gA"^ = 2.75, and a bare quark mass m = 5 MeV. With these 
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Figure 1: Various quantities obtained with the Lorentz invariant interaction = gu = 2.75/A^ 
as functions of the quark chemical potential /x. Left panel: (solid), Ai (dashed), /x^ (dashed- 
dotted), /ife (dotted). Right panel: — Mi, (solid), — A2 (dashed). 



parameter values we obtain a vacuum constituent quark mass Mr = Mi, = 407.7 MeV. 
This corresponds to a quark condensate = — 2(245. 7MeV)^, while 08, P, Ps, and 62 
vanish in vacuum. 

When we increase the quark chemical potential nothing happens up to a critical value 
Pcrit = 403.3 MeV. At this point a first-order phase transition takes place and all expecta- 
tion values Eqs. (|l|) to (|^) receive non-vanishing values. This can be inferred from Fig. |l] 
where various quantities are displayed as functions of /i. In the left panel we show the con- 
stituent quark mass M,., the diquark gap Ai, and the effective chemical potentials fir and 
/ift. In the right panel the mass difference Mr — M^ and the diquark gap — A2 are plotted. 
At /i = ficrit the constituent quark masses drop by more than 300 MeV and are no longer 
identical. The difference, however, is small, Mr = 95.7 MeV and Mb = 95.0 MeV. With 
increasing chemical potential, both, the masses and their difference, are further decreased. 
In the diquark channel we find Ai = 120.0 MeV at = Pcrit- Similar to what has been 
found in Ref. 0, the second gap parameter has the opposite sign and is more than one 
order of magnitude smaller, A2 = -8.4 MeV. Like the constituent masses, it decreases 
with increasing n, whereas Ai is slightly growing in the regime shown in Fig. |^. 

Below the phase transition, the effective chemical potentials fir and fif, are equal to the 
external chemical potential fi. At the phase transition point, fir and fib drop by 67 MeV 
and 51 MeV, respectively and then grow again as functions of fi. The corresponding 
number densities of red and blue quarks are shown in Fig. |[ At /z = pcrit, Pr jumps 
from zero to 0.42 fm~^, about 2.5 times nuclear matter density, whereas the density of 
blue quarks reaches only twice nuclear matter density at this point, pb = 0.34 fm~'^. Both 
densities grow of course with increasing chemical potential, but their difference remains 
nearly constant. 

The unequal densities of red and blue quarks in this state, although anticipated, clearly 
deserve further discussion. First of all we should point out that pr — pb ^ is not a result 
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Figure 2: Left: Number densities of red quarks (solid) and blue quarks (dashed) as functions of 
the quark chemical potential fi. Right: Energy per quark as function of the total quark number 
density for a color superconducting system with equal densities of gapped and ungapped colors 
(dashed) and with unequal densities as given in the left panel (solid). 



of our self-consistent treatment, but is already found if we only take into account the 
effect of a non- vanishing Ai. Nevertheless, it is usually not discussed in the literature. 

Since our model Lagrangian Chge is symmetric under global color SU{3) transforma- 
tions, the total number of quarks of each color is a conserved quantity. As in BCS theory, 
quark number conservation is violated by the diquark condensates, but the average num- 
bers are still conserved and given by the densities pi integrated over the volume. One 
could therefore ask what happens if we start with a large but finite system with equal 
numbers of red, green and blue quarks at low densities and then compress it until a color 
anti-triplet diquark condensate is formed. According to the above results, in this phase 
there are not equally many quarks of each color, but the number of those quarks which do 
not participate in the condensate (in the above case the blue ones) is smaller. Obviously, 
we cannot get a single phase with these properties if we start from a system with equal 
densities. A possible scenario could be that several domains emerge in which the symme- 
try is broken into different directions, such that the total number of red, green and blue 
quarks remains unchanged. Still one could worry about strong forces inside these domains 
arising from the non-vanishing net color charge of the domain. However, we should keep 
in mind that by construction our quarks are non-interacting quasiparticles, i.e., at least 
on the mean field level there are no forces left. 

Alternatively, we can construct a superconducting state with equal densities for the 
gapped and ungapped quarks. To that end we have to introduce different external chem- 
ical potentials for different red and blue quarks, or, equivalently, an additional external 
chemical potential /ig^*. Then the second equation in Eq. (|Tip becomes 

/is = /i^'* + 2(7(^^8 , (36) 
and in Eq. (|1^) we should replace /ig by {ps ~ I^T^Y- With this additional external pa- 
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rameter we can enforce the densities of all colors to be equal, even in the superconducting 
state. Obviously this is the case if ps = f^T^- 

However, at least within our mean field approximation such a state would be energet- 
ically less favored than a state with the same total density, but = 0. This is shown 
in the right panel of Fig. 0. The energy density e of the system is given by 

e{T = 0, p, ps) = n{T = 0, p, ps) - ^](0, 0, 0) + pp + pf p8 . (37) 

The second term on the r.h.s. has been introduced to normalize the energy density of 
the (non-trivial) vacuum to zero. In the right panel of Fig. ^ the energy per quark, e/p, 
is displayed as a function of the total quark number density p. The solid line is the 
result for pg"^* = 0, i.e., it corresponds to the unequal red and blue quark densities as 
shown in the left panel. The dashed line corresponds to the result for equal red and 
blue quark densities. As one can see in the figure, the energy of this solution is always 
higher than the energy of the solution with unequal densities. This means, according to 
this result, a large homogeneous system of equally many red, green and blue quarks is 
unstable against decay into several domains in which the density of the gapped quarks is 
larger than the density of the remaining ungapped quark. On the other hand, the energy 
difference is not very large (less than 1 MeV per quark). Therefore it cannot be excluded 
that the homogeneous solution with equal densities would be favored, once effects beyond 
the mean- field approximation are taken into account. 



4.2 Lorentz non-invariant interaction 

In the previous section we found rather small values for the gap parameter A2. This seems 
to justify the common practice to neglect this condensate completely. In order to find out 
whether the smallness of A2 is a general feature, we now want to abandon the unnecessary 
restriction to a Lorentz invariant interaction, Eq. (p5|), and to allow for unequal values 
oi gE and gM- Looking at Eq. (|3l), we see that the coupling constants gf'\ gi^'^ and hi 



depend on the sum gE + SgM, whereas gv^^ and /i2 depend on the difference gs — SgM- 
gi is somewhat exceptional because it is mostly given by the "direct" interaction term, 
which is not present in the other channels. The "exchange" term is again proportional to 
gE — SgM, but suppressed by a factor 1/24. 

In principle there are two ways to obtain larger values of A2. The more dramatic one 
is to assume that gE is larger than 3gM- Then the coupling constants glj'^ and /i2 are 
positive, i.e., the interaction becomes attractive in these channels. Indeed, as we will see 



in Sec. |4.3| , an attractive /i2 can yield solutions with large values of A2. However, these 
solutions might not be stable, similar to what we found for the schematic example in the 
end of Sec. 0. Unfortunately, a detailed study of the structure of solutions of the coupled 
gap equations involving four attractive channels is very difficult and beyond the scope of 
our paper. 

The second possibility is to make h2 more repulsive. In order to have a similar starting 
point as in the previous section, we keep the sum gE+'igu fixed, but increase the difference 
Sfi'M — gE- To see a relatively large effect we take gE = — 1.75/A^ and gM = 4.25/A^, i.e., 
Sfi'M — gE = 14.5/A^, almost three times as large as in the previous section. As before, we 
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Figure 3: Various quantities obtained with the Lorentz non-invariant interaction, = 
— 1.75/A^ and qm = 4.25/A^. Left panel: Mr (solid), Ai (dashed), — A2 dotted, and fir (dashed- 
dotted) as functions of the quark chemical potential Right panel: Dispersion law E-{p) at 
/i = 450 MeV. The dashed line was calculated neglecting A2 in the gap equations, whereas the 
solid line corresponds to the exact solution. 



take m = 5 MeV and a cutoff A = 600 MeV. Our results are displayed in Fig. ^ In the left 
panel we show the behavior of Mr, Ai, — A2, and fir as functions of fi. The most striking 
difference to our previous example (Fig. |ID is the fact that we now find a smooth crossover 
instead of a first-order phase transition. In the chiral limit the phase transition becomes 
second order. A similar effect is known from studies of the chiral phase transition within 
the NJL model without diquark condensation. There it was also found that a first-order 
phase transition becomes second order (or a smooth crossover) if the coupling strength in 
the (repulsive) vector channel exceeds a certain value [jl6|. Although the present model 



is rather unrealistic below n ~ 450 MeV, where it predicts a color-superconducting quark 
gas of low density, it is nevertheless a counter example to the common (and mostly model 
based) belief that the chiral phase transition at zero temperature and large n should be 
first order (see e.g. Refs. |Q, p^). 

Despite the fact that |/i2| is larger than hi, the absolute value of A2 remains always 
smaller than that of Ai. However, unlike in the previous section, A2 can no longer 
be generally neglected. For instance, at /i = 450 MeV, we find Mr = 280.1 MeV, 
Hr = 382.3 MeV, Ai = 83.2 MeV, and A2 = -37.8 MeV. In terms of Eqs. (^ and 
(13) this corresponds to /ie// = 384.1 MeV, M^// = 286.9 MeV, and |Ae//| = 55.2 MeV. 
The resulting dispersion law E-{p) is shown in the right panel of Fig. ^ (solid line). At 

\p\ = ^JfJ^lff - Ml^j = 255.4 MeV it has a minimum with E_ = |Ae//| = 55.2 MeV. On 
the other hand, if we neglect A2 in the gap equation, we get M^ff = Mr = 292.5 MeV, 
fieff = fJ^r = 383.3 MeV, and Ae// = Ai = 102.2 MeV. Consequently, the minimum value 
of E_ is now 102.2 MeV, almost twice as much as without neglecting A2. The corre- 
sponding dispersion law is indicated by the dashed line in the right panel of Fig. 0. As 
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one can see, the entire function E^{p) is shifted to higher energies as compared with the 
sohd curve and the minimum is more shallow. 



4.3 Gapless color superconductors 

Finally, we would like to come back to the question of possible gapless color superconduc- 
tors. Obviously, none of our numerical examples presented so far came close to condition 
(H). For instance, if we take Mr = 280.1 MeV, /i^ = 382.3 MeV, and Ai = 83.2 MeV, as 
found in the previous section at fj, = 450 MeV, one would need A2 — -113.6 MeV, about 



three times as large as the actual value. The situation was even worse in Sec. |4.1| where 
the discrepancy was about a factor 50 at /i = ficrit and became larger with increasing 
chemical potential. In fact, none of our numerical examples fulfilled IA2I > |Ai|, which 
we identified as a necessary condition for gapless color superconducting states. 

To get some insight, how an interaction could look like which yields such a state, we 
can invert the problem and employ the gap equations to calculate the effective coupling 
constants which are consistent with a given set of gap parameters. For instance, Eq. (|2^) 
is obviously fulfilled if we choose = Ai = 100 MeV and Hr = -A2 = 350 MeV. For 
simplicity we might assume Mj, = Mr and /i^ = fir- Except of A2 this is within the 
typical range of these quantities in the earlier examples. If we now take m = 5 MeV and 
a cutoff A = 600 MeV, as before, and /i = 450 MeV, the gap equations yield gf'^A'^ = 
3.36, gi^^A'^ = -1.41, /iiA^ = 6.80, /laA^ = 6.18, and gi^^A^ = gi^^A^ = 0. Here the 
essential difference to our earlier examples is the need of an attractive interaction in the 
/i2 channel. Furthermore, the interaction is relatively strong in both diquark channels. 
However, for these parameters there is another solution with Mr = Mb = 58.1 MeV, 

= ^f, = 362.8 MeV, Ai = 966.6 MeV, and A2 = 16.1 MeV. In order to decide which of 
the two solutions is the energetically favored one we have to evaluate the thermodynamic 
potential. It turns out that for the gapless solution the value of Q is about 900 MeV/fm^ 
higher than for the other solution. Hence, similar to what we found in the schematic 
example of Sec. |^, the gapless state does not correspond to a stable solution. In fact, we 
did not succeed to construct a stable color-superconducting solution. This might indicate 
that such a state does not exist, although a rigorous proof is still missing. 



5 Conclusions 

We analyzed the ground state properties of an isotropic two-flavor color-superconductor 
at finite quark chemical potential within an NJL-type model. In general, a self-consistent 
treatment of this problem requires to consider several condensates which are usually ne- 
glected: The breakdown of Lorentz invariance in dense systems implies the possible ex- 
istence of Lorentz non-invariant diquark condensates, while the spontaneous breaking of 
color SU{3) in a color superconductor naturally leads to the existence of SU{3) non- 
invariant quark-antiquark condensates. We found that at least six different expectation 
values (two diquark condensates and four quark-antiquark condensates, Eqs. to (^)) 
have to be taken into account in a self-consistent calculation and we derived a set of six 
coupled gap equations for these expectation values. 
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The actual importance of the various condensates depends of course on the interaction. 
Since at present not very much is known about the effective interaction which describes 
the deconfined phase at moderate densities, this leaves room for surprises and possible 
new phenomena. For instance, in the numerical example discussed in Sec. [4.1| the chiral 
phase transition at finite fi and zero temperature was of first order, in agreement with 
the general expectation. However, in Sec. |4.2| , where a Lagrangian with strong repulsive 
interactions was chosen, we found a smooth crossover. Similarly, the effect of the Lorentz 
non- invariant diquark condensate 62, which is usually neglected, was indeed found to be 
small in Sec. [4.1| , but relatively large in Sec. f4.2| . There it caused a reduction of the effective 
gap parameter |Ae//|, (the minimum energy of a quasiquark excitation) by almost 50%. 
In general, the numerical values of the condensates depend also on how the divergent 
integrals are handled. In this article we used a sharp cutoff. A comparison with other 
approaches is yet to be done. 

We also discussed the possibility of "gapless color superconductors", i.e., states with 
non- vanishing diquark condensates, but |Aejj| = 0. We have seen that such states can 
in principle exist, provided the gaps of the various condensates are related to each other 
in a certain way (see Eq. (0)). In practice, however, all gapless superconducting states 
we constructed turned out to be unstable solutions of the gap equations. A similar 
observation was made in Ref. [|15| for gapless states in the color-flavor locked phase. This 
suggests that gapless color superconductors might be in general unstable. 

As a consequence of the spontaneous color SU (3) breaking through the diquark con- 
densates, there are also SU{3) non-invariant quark-antiquark condensates in a color su- 
perconductor. For instance, the familiar scalar condensate (ipip) is in general different for 
quarks which participate in the diquark condensate ( "red" and "green" ) and those which 
do not ("blue"). This also leads to different constituent quark masses for gapped and 
ungapped quarks. However, at least in our numerical examples this difference turned out 
to be quite small. On the other hand we found that (for equal chemical potentials) the 
density of the gapped quarks was considerably larger than the density of the ungapped 
quarks. Since the total number of red, green and blue quarks should be equal in a finite 
system, this could lead to the emergence of domains in which color SU{3) is broken into 
different directions. 
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